Effects of Post-Newtonian Spin Alignment on the Distribution of Black-Hole Recoils 
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Recent numerical relativity simulations have shown that the final black hole produced in a binary 
merger can recoil with a velocity as large as 5, 000 km/s. Because of enhanced gravitational-wave 
emission in the so-called "hang-up" configurations, this maximum recoil occurs when the black-hole 
spins are partially aligned with the orbital angular momentum. We revisit our previous statistical 
analysis of post-Newtonian evolutions of black-hole binaries in the light of these new findings. We 
demonstrate that despite these new configurations with enhanced recoil velocities, spin alignment 
during the post-Newtonian stage of the inspiral will still significantly suppress (or enhance) kick 
magnitudes when the initial spin of the more massive black hole is more (or less) closely aligned 
with the orbital angular momentum than that of the smaller hole. We present a preliminary study 
of how this post-Newtonian spin alignment affects the ejection probabilities of supermassive black 
holes from their host galaxies with astrophysically motivated mass ratio and initial spin distributions. 
We find that spin alignment suppresses (enhances) ejection probabilities by ~ 40% (20%) for an 
observationally motivated mass-dependent galactic escape velocity, and by an even greater amount 
for a constant escape velocity of 1,000 km/s. Kick suppression is thus at least a factor two more 
efficient than enhancement. 

PACS numbers: 04.25.dg, 04.70.Bw, 04.30.-w 



I. INTRODUCTION 

Supermassive black holes (SBHs) reside at the centers 
of most large galaxies. The masses of these SBHs are 
tightly correlated with the luminosity [1[, mass [2|, and 
velocity dispersions 0] of the spheroidal components of 
their host galaxies. Large galaxies form hierarchically 
from the merger of smaller galaxies, and SBH mergers 
are expected to accompany the mergers of their hosts. 
The final stages of these SBH mergers are driven by the 
emission of copious amounts of gravitational radiation. 
Conservation of linear momentum implies that the final 
black holes produced in SBH mergers must recoil with 
linear momentum equal in magnitude and opposite in 
direction to that of the anisotropically emitted gravita- 
tional radiation. Early estimates of these gravitational 
recoils or "kicks" using post-Newtonian (PN) techniques 
[J] and black-hole perturbation theory [5[ suggested that 
they would not exceed several hundred km/s in mag- 
nitude. More recently, progress in numerical relativity 
(NR) p4al has allowed relativists to simulate the merg- 
ers of highly spinning, comparable-mass black holes. For 
non-spinning binaries, ensuing studies identified a maxi- 
mum recoil of 175 km/s [9|. Simulations of spinning black 
holes, however, resulted in one of the greatest surprises 
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of numerical relativity; for equal-mass binaries with op- 
posite spins in the orbital plan e gravitational recoils can 
approach 4,000 km/s [ifl, ^^, greater than the escape 
velocity of even the most massive galaxies. This theo- 
retical result seems at first difficult to reconcile with ob- 
servations indicating that almost all large galaxies host 
SBHs. 

One solution to this problem is to align the SBH spins 
before merger into configurations that lead to small re- 
coils. Gravitational radiation extracts energy and angu- 
lar momentum from the orbit of the binary SBHs, causing 
them to merge on a timescale <gr cx: r'^, where r is the 
binary separation. This timescale becomes longer than 
the age of the universe at separations r ~ 1 pc, implying 
that some mechanism other than gravitational radiation 
is required to escort the SBHs to this separation (the 
"final-parsec problem" ) [12| . One such mechanism is the 
transfer of orbital angular momentum from the SBH bi- 
nary to surrounding gas. This gas forms a circumbinary 
disk about the SBHs, and can be transferred through the 
gap onto accretion disks about the individual SBHs [ij] . 
If the angular momentum of the circumbinary disk is 
misaligned with the SBH spins, the Lense-Thirring effect 
causes inclined annuli in the individual accretion disks to 
differentially precess about the SBH spins. Bardeen and 
Petterson [lj| showed that viscous dissipation causes this 
differentially precessing gas to settle into the equatorial 
planes of the SBHs. On longer timescales, these warped 
accretion disks torque the SBH spins into alignment with 
the orbital angular momentum of the gas at large radii, 
presumably that of the circumbinary disk from which 



both accretion disks arc being fed [l5|. Bogdanovic et 
al. [l^ suggested that this alignment could reduce grav- 
itational recoils to less than 200 km/s, in contrast to the 
^ 4, 000 km/s recoils expected for SBHs in the "super- 
kick" configuration (s pin s in opposite directions lying in 
the equatorial plane) [10, llll ■ 

Dotti ct al. [131 tested this suggestion by perform- 
ing a series of N-body smoothed particle hydrodynamics 
(SPH) simulations of two 4 x 10^ Af© SBHs inspiraling 
due to dynamical friction exerted by a IO^Mq circumnu- 
clcar disk in their orbital plane. One of the SBHs began 
at the center of the circumnuclear disk, while the second 
SBH spiraled inwards from an initial separation of 50 pc 
to a final separation of ~ 10 pc. Gas particles within the 
Bondi-Hoyle radii [l3l of the SBHs were accreted, and 
assumed to fuel warped accretion disks as described in 
Perego et al. [il|. On a timescale of < 1-2 Myr, the SBH 
spins became aligned to within 10° (30°) of their orbital 
angular momentum for a cold (hot) circumnuclear disk 
[l9[ . If the partially aligned spin configurations found in 
these simulations were preserved as the SBHs inspiralcd 
and merged, the median recoils predicted by NR simu- 
lations would be < 100 km/s for dimensionless spins Xi 
[i = 1,2) with an initially uniform distribution of mag- 
nitudes Xi G [0, !]• 

The assumption that SBH spin distribu- 
tions remain unchanged between r ~ 10 pc ~ 
lO^Rs (Mbr/IO^Mq)-^ (where the SPH simulations 
end) and r ~ 5 Rs (where the NR simulations begin), 
needs further examination. Here Rs = 2GMbh/c^ is the 
Schwarzschild radius of a SBH of mass Mbh- Because 
of the steep scaling tcR oc r^ mentioned previously, 
the SBHs decouple from their circumbinary disk and 
spiral inwards purely under the influence of gravitational 
radiation at a binary separation rdec — 10^ — 10'^ M 
jl2| . where M = mi + m2 is the total mass of the 
binary. Schnittman [2^ integrated PN equations of 
motion and spin precession from an initial binary 
separation r^ = lOOOM to a final separation r/ = lOM 
(in geometrical units where G = c = I) and discovered 
that partially aligned SBH spin distributions during this 
portion of the inspiral are strongly influenced by the 
presence of spin-orbit resonances. 

These resonances are special spin configurations in 
which both SBH spins and the orbital angular momen- 
tum jointly precess at the same frequency in a common 
two-dimensional plane. For binary SBHs with mass ratio 
q = m2/mi < 1 and dimensionless spins Xi: ^t each sep- 
aration r there exist two one-parameter families of spin- 
orbit resonances: one with A0 = 0° and the other with 
A0 = 180°, where A(j) is the angle between the compo- 
nents of Xi s-nd X2 perpendicular to the orbital angular 
momentum L. The first of these families (A0 = 0°) has 
di < 02, where dj = arccos(Xj • L), while the second fam- 
ily {A(t> = 180°) has 6*1 > 6*2. As SBHs inspiral due to 
the emission of gravitational radiation, they can become 
captured into nearby spin-orbit resonances, substantially 
altering the SBH spin distributions between r^ and rf. 



In previous work, we examined the influence of this 
PN portion of the inspiral on the expected distribution of 
SBH final spins [2l| and recoils [Hj. We found that if the 
spin of the more massive SBH was more closely aligned 
with the orbital angular momentum {9i < O2), the SBHs 
were more likely to be captured into the A0 — 0° fam- 
ily of resonances. This alignment of the SBH spins prior 
to merger causes them to add constructively, enhancing 
the spin of the final black hole. Symmetry requirements 
imply that gravitational recoils are suppressed when the 
SB H sp ins are aligned with each other prior to merger 
[23I |24| , so this PN alignment also reduces the predicted 
distribution of recoils. The opposite is true when 6*1 > 6*2; 
SBHs are preferentially captured into the A<j> ~ 180° 
family of resonances, reducing the final spins and enhanc- 
ing the gravitational recoils. 

Recen t eq ual-mass NR simulations by Lousto and Zlo- 
chower ([2^, [2g; henceforth LZ) have qualitatively and 
quantitatively changed the predicted dependence of grav- 
itational recoils on the SBH spins Xi- The maximum pos- 
sible kick, extrapolated to maximal initial spins, is now 
~ 4, 900 km/s, and this kick occurs not in the previously 
described "superkick" configuration but for spins with 
9i ~ 50°. In light of these new findings, it is worth inves- 
tigating whether PN spin alignment still has as dramatic 
an effect on the expected recoil distribution as we found 
previously. This investigation is the subject of this paper. 
In Section |TT] we revisit our analysis of how PN evolution 
affects binaries that are evolved from an initial separation 
of Vi = lOOOM down to a final separation of rf = lOM. 
We clarify the reason why spin-orbit resonances are also 
effective at suppressing kicks consistent with the new LZ 
results. In Section IIIII we set up a more extensive set 
of Monte Carlo simulations of SBH inspirals. We focus 
on a specific choice for the initial distribution of SBH 
parameters {9, Xi,X2}i ^^'^ justify this choice based on 
astrophysical considerations. In Section IIVI we present 
the final distributions of recoil velocities resulting from 
our Monte Carlo simulations, and our predictions for the 
fraction of remnant SBHs that are ejected from their host 
galaxies. Wc summarize these results and discuss the 
sources of uncertainty in our analysis in Section [V] 



II. POST-NEWTONIAN SPIN ALIGNMENT 

The recent LZ simulations [2^,[26| have two potentially 
important consequences from an astrophysical point of 
view. First of all, the maximum possible "hang-up" kick 
magnitude is about 25% larger than predicted by the 
ordinary, "no hang-up" superkick formula: 4,900 km/s 
rather than 3,750 km/s. Secondly, the largest kick occurs 
for a configuration in which the spins are partially aligned 
with the orbital angular momentum. Because of the gas- 
dynamical alignment arguments summarized above, this 
configuration seems more likely than the superkick con- 
figuration. Wc would therefore expect an enhancement 
in the statistical likelihood of gravitational recoils eject- 



ing remnant SBHs from galaxies as reported by LZ. We 
confirm this enhancement, but show that even in hang- 
up kick scenarios spin-orbit resonances can significantly 
change the likelihood of large recoil velocities as the SBHs 
inspiral from sub-parsec scales to about lOM. 

The largest "superkick" recoil velocities are obtained 
for nearly equal-mass SBHs. Furthermore, the LZ kick 
formula is most reliable for mass ratios close to unity, 
since all 48 of the new NR simulations were carried out for 
q = I [2al ■ We therefore begin our study of the likelihood 
of very large kicks by considering comparable-mass bina- 
ries. Resonance effects are not present for q = 1, but as- 
trophysical binaries are not expected to be precisely equal 
in mass. For this reason we have performed Monte Carlo 
PN evolutions of 900 black hole binaries with mass ratio 
close to but not exactly equal to one (we choose q = 9/11 
to facilitate comparisons with Schnittman [20[ and our 
own earlier work [2l|, ^^ ) . We considered binaries where 
the two black holes have the same dimensionless spin 
magnitude Xi = X2 = Xi and selected five equally-spaced 
values of x in the range [0.5, 1.0]. For each of the five 
values of x we further selected six equally-spaced initial 
values of 6i in the range [5°, 30°] and six values in the 
range [150°, 175°], for a total of sixty Monte Carlo simu- 
lations. Each simulation is started by assuming uniform 
distributions in A(f> £ [0 ,2tt] and cos 02 G [^Ijl]- We 
begin at an initial separation r^ = lOOOM and use the 
PN equations of motion to evolve the binaries down to 
a final separation r/ = lOM (cf. Section II of [2l| for 
details of the PN evolution). 

To quantify the effect of PN spin alignment, wc com- 
pute the recoil velocity distribution that would result 
from applying kick formulae to our binaries before the 
PN evolution (at r^ = lOOOM) and close to merger (at 
Tf = lOAf). We calculate the component W|| of the recoil 
velocity par allel to L using both the new LZ "hang-up" 
formula [251, 1261 and the older kick formula from Cam- 
panelli et al. [ll|. We calculate the components Vm, 
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v± perpendicular to L from the fitting formula of 
cf. Eq. (1) in [221 ■ We then compute the probability that 
the total recoil velocity is greater than v^j = 1, 000 km/s 
for each of our sixty Monte Carlo simulations (twelve 
values of 9i and five values of x)- This value of Woj is 
comparable to the escape velocity from a giant elliptical 
galaxy [23| ; we will consider more realistic escape veloci- 
ties that depend on the host galaxy's mass in Section HVl 
The results arc plotted in Figure [TJ 

Figure [1] confirms LZ's conclusion that the "hang- 
up" effect increases ejection probabilities if the spin of 
the heavier companion is partially aligned with the or- 
bital angular momentum (^i < 30°) whereas the ejection 
probabilities are reduced for partial antialignment of xi 
(6*1 > 150°); cf. their Fig. 2. By comparing the dashed 
lines in the upper left and right panels of Figure [3 we 
see that, if we ignore the effect of resonances (as the au- 
thors of the LZ papers did), ejection probabilities for the 
aligned cases 9i < 30° are only mildly dependent on 6i 
and roughly double because of the hang-up kick effect. 
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FIG. L Probability of having a recoil velocity v > 1, 000 km/s 
as a function of dimensionless spin magnitude Xi — X2 = X 
for binaries with mass ratio q = 9/11. The initial value of 
9i (the alignment angle of the larger SBH) is indicated by 
the different symbols, while the spin X2 '^^ the smaller SBH 
is initially isotropically distributed. Left panels: prediction 
according to the LZ papers [251. 1261]: right panels: prediction 
without this newly discovered hang-up effect. In all plots the 
dashed lines correspond to the initial distributions at n = 
lOOOAf , while the solid lines give the distributions at r/ = 
lOAf , after PN spin alignment has occurred. 



These ejection probabilities can be larger than 50% for 
spins X ^ 0.8 and q — 9/11. On the other hand, a 
comparison of dashed and solid lines in each panel shows 
that ejection probabilities are dramatically reduced in the 
aligned case if we apply the kick formula (as we should) 
close to merger, rather than at large separations. The 
relevance of resonant effects depends on 9i : if the initial 
angle between the larger black hole and the orbital an- 
gular momentum 9i < 10°, recoil velocities larger than 
Doj ^ 1,000 km/s would almost never occur. Even for 
more modest alignments {9i < 30°), such large kicks 
only occur for large values of x- Recoils are even smaller 
for q < 9/11, as all components of the kick are propor- 
tional to q'^/il + q)^ (see Eq. (11) of (Hi). Note that this 
kick suppression only applies when on average 9i < 92', 
for initial values of 9i > 6*2 the recoil is enhanced, not 
suppressed; as demonstrated by the lower two panels of 
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FIG. 2. Left panel: histograms of the recoil velocity v for maximally spinning (|Xil = 1) BH mergers with mass ratio q = 9/11. 
Dotted lines show the recoils expected if the BHs merge with the parameters specified at initial separation Vi — lOOOAf . Solid 
curves show the expected kicks if the BHs precess as described in Section 2 of [22 from r^ to r/ — WM prior to merger. Right 
panel: same as in the left panel, except that here the solid histograms show the velocity distribution obtained when Acp is reset 
to a flat distribution A(j) £ [0, 180°] after the binaries reach r/. 



Fig.[TJ where 9i > 150°. Interestingly, the kick enhance- 
ment appears to saturate: while ejection probabilities ap- 
proach zero as we fix 6i at ever smaller values, they are 
essentially independent of 150° < Oi < 175° if we apply 
the kick formula close to merger. Figure [T] also indicates 
that kick suppression via alignment {6i ~ 0°) is more ef- 
ficient than enhancement via antialignment (^i ^ 180°) . 
We shall quantify this asymmetry in more detail below. 

It may seem surprising that PN spin alignment so dra- 
matically suppresses the recoil velocities predicted by the 
new LZ kick formula, since many of the spin-orbit res- 
onances depicted in Figure 1 of [21| have values of 9i 
comparable to that of the new "maximum recoil" con- 
figuration in the presence of the hang-up effect. We can 
understand the effectiveness of PN kick suppression by 
recognizing that it results from the alignment of the per- 
pendicular components of the SBH spins with each other 
(Acf) — !■ 0°), not from the spins aligning with the orbital 
angular momentum (0^ — >■ 0, as Bogdanovic et al. [1^1 ar- 
gued would occur due to the torque exerted by a warped 
circumbinary disk). 

To illustrate this point, we consider maximally spin- 
ning (xi = X2 = 1) binaries with q = 9/11 in Figure [2j 
Dotted histograms refer to recoil distributions computed 
before the PN evolutions (at r^ = lOOOM), and solid 
histograms refer to the corresponding distributions com- 
puted after evolving the binaries down to rf = lOil/. 
Let us first focus on the left panel. The black curves in 
the upper histogram show that black holes with isotropic 
spin distributions (flat distributions in cos 9i , cos 02 , and 
A(/)) maintain these isotropic distributions as they in- 
spiral down to rf = IQAI. The thick blue (thin red) 
curves in the upper histogram corresponds to the subset 
of this isotropic distribution with the 30% lowest (high- 



est) initial values of 9i. Careful comparison of the solid 
and dotted curves reveals that recoils are suppressed (en- 
hanced) for low (high) initial values of 6i . This tendency 
is seen much more clearly in the middle and lower panels, 
where 9i has the indicated initial value while X2 retains 
its isotropic initial distribution. If Si is partially aligned 
with L, as in the middle histograms, then ^12 (the angle 
between the two spins) and A4> will be strongly peaked 
around 0° at the end of the evolution and kicks will be re- 
duced. If instead Si is partially antialigned with L (as in 
the lower histograms), then 9i2 and A(f> will be strongly 
peaked around 180° at the end of the evolution and the 
kicks will be enhanced (cf. the discussion in [22[)- The 
crucial element for producing alignment here is the fact 
that A(f> -^ 0°. This is illustrated in the right-hand pan- 
els of Figure [U which are the same as the corresponding 
left-hand panels, except that A(j> has been reset to a flat 
distribution A(j) e [0, 180°] after the PN inspiral (but be- 
fore the recoil velocities are computed) . This nearly elim- 
inates the kick suppression (for initially aligned binaries) 
or enhancement (for initially anti-aligned binaries). 



III. ASTROPHYSICAL PARAMETERS 

We demonstrated in the previous section that the new 
fitting formulae for gravitational recoils provided in the 
LZ papers [25|, |26[ did not alter our earlier conclusion 
[23 ] that PN spin alignment can dramatically affect the 
distribution of recoil velocities for specific initial values 
of q and 61 and an initially isotropic distribution for X2- 
In this section, we describe a more astrophysically mo- 
tivated choice for these initial parameters. The effect of 
PN spin alignment on this new distribution will be pre- 



sented in Section IrVI 



by stellar mass; we use the functions 



A. Mass Ratio 

The origin of SBHs is poorly understood theoretically 
and poorly constrained observationally. Lynden-Bell [28| 
recognized that the central objects of mass 10^ — 10^ Mq 
believed to power quasars would quickly collapse into 
SBHs. Haehnelt and Rees [23 explained the quasar lumi- 
nosity function by assuming that such massive SBH seeds 
promptly formed at the centers of dark-matter (DM) ha- 
los whose formation rates could be predicted using the 
Press-Schechter formalism [301 ■ More recent theoretical 
work has called into question whether the high-redshift 
seeds of SBHs were truly this massive. While some main- 
tain that massive central concentrations of gas at the 
centers of pre-galactic disks can directly collapse into 
~ 10^ Mq SBH seeds [3l|, others argue that metal-free 
gas will fragment into Population HI stars with charac- 
teristic masses > IOOMq [32|. Volonteri, Haardt, and 
Madau [Sj] , also using the Press-Schechter formalism to 
estimate halo mass functions and merger rates, showed 
that 150 Mq seeds occupying 3.5 — 4cr overdensities at 
redshift z — 20 could grow into SBHs that would repro- 
duce both the quasar luminosity function at 1 < z < 5 
and the observed local A/bh — o" relation [^ between SBH 
mass and galaxy velocity dispersion. A future space- 
based interferometer such as LISA or the European New 
Gravitational Wave Observatory (NGO/eLISA) will be 
able to distinguish between these two different scenarios 
for SBH formation by observing the grav itational waves 
emitted during high-redshift mergers (34| - [37| . 

Given this uncertainty in the high-redshift population 
of SBH seeds, in this paper we restrict our attention to 
SBH mergers in the low-redshift (z < 1) universe. Stew- 
art et al. [331 used a high-resolution N-body simulation 
to calculate the merger rates of DM halos in the red- 
shift range < z < 4. They fit their results for the rate 
at which a larger galaxy of mass Mg merges with smaller 
galaxies with masses between nig and Mg to the function 



dN, 



incrg 



dt 



{Mg,mg)^At{Mg)F{mg/Mg). 



(1) 



where the normalization At{Mg) is binned by mass, and 
the mass-ratio dependence is given by 



(2) 



i^K/M.) = hrf i-TT 



Mn 



Ma 



with the indices c and d determined from the simulation. 
This same functional form can be used for the merger 
rates of both DM halos of mass Mh and galaxies with 
stellar mass Af*, although the fitted values of the pa- 
rameters will be different. Table 1 of Stewart et al. [S^ 
provides numerical estimates of these parameters binned 



At(M^) = 0.0098 Gyr" 



M, 



c = 0.329 



A4 



IQi" Mq 
0.158 



0.736 



d= 1.1-0.2Z, 



, (3a) 

(3b) 
(3c) 



which approximate the values listed in this table. Note 
that the total rate of mergers diverges in the limit mg — > 
for c > 0, although the rate at which the galactic 
mass increases remains finite. The Press-Schechter for- 
malism also predicts that the number of DM halos and 
thus the number of mergers diverges as their mass goes 
to zero. Volonteri, Haardt, and Madau [SJl address this 
issue in their merger trees by only keeping track of halos 
with masses above an effective mass resolution Aires that 
evolves with redshift, but always remains a small fraction 
of the largest progenitor mass. Because we are only inter- 
ested in SBH mergers at low redshift, we neglect galax- 
ies with Af* < 10^ Af©. This corresponds to neglecting 
SBHs with nii < 10^ Mq, since SBHs are about a fraction 
of 10^"^ of the stellar mass of their host spheroids [39| . 
SBHs less massive than this value are poorly constrained 
observationally, and even if present in a halo may have 
difficulty reaching the galactic center, since the dynam- 
ical friction time scales as the inverse of the SBH mass. 
We assume that SBH mergers occur promptly following 
the merger of their host galaxies. 

Eq. ID) gives the rate at which individual galaxies of 
stellar mass Af* merge with smaller galaxies; to deter- 
mine the total rate of mergers in a cosmological volume 
also requires an estimate of the galactic stellar mass func- 
tion dNg/dlogM^,. We use the mass function determined 
from the z ~ 0.37 sample of the COSMOS survey ^; 
this mass function, as shown in the left panel of Figure 14 
of [iol, is quite similar to earlier mass functions derived 
from the Sloan Digital Sky Survey filj - lisj . The total rate 
per unit cosmological volume i?(Af*) at which galaxies 
with stellar masses between 10^ Mq and Af* merge is 
given by 



i?(A/,) 



'°'''' ^^^%Ii(Af^, io^Mq) dlogMl . 
dlogA/. dt ^ *' oj a * 

(4) 
We use this rate function as part of a Monte Carlo ap- 
proach to generate a representative sample of SBH merg- 
ers. We choose a random number ri from a flat distri- 
bution in the interval < ri < 1, then find the stellar 
mass M^. for which R{M^) = 7'ii?(oo). Values of ri ~ 
correspond to stellar masses A/* ~ 10^ Mq, while values 
ri ~ 1 correspond to stellar masses Af* ~ IO^^Mq above 
which the galactic luminosity function is exponentially 
suppressed. We next choose a second random number 
r2, also from a flat distribution < r2 < 1, and find the 
value of m* for which F{m^/M^) ^ r2F {10^ Mq/M^). 
Values of r2 ~ correspond to nearly equal-mass mergers 
(jTi* ^ Af*), while values r2 — I correspond to mergers 



with the smaUest galaxies to host SBHs (m* ~ 10^ Mq). 
The SBH masses are then determined from the masses 
of their host galaxies, nii — 10~^M* and 7712 = 10~'^to*. 
Histograms of our SBH mass ratio distribution are given 
by the solid black curves in the upper panels of Fig. [3] 

Our assumption that the SBH mass is 10""^ that of its 
host galaxy's stellar mass is strictly true only for elliptical 
galaxies, since observations suggest that SBH masses are 
uncorrelated with galactic disks and pseudobulges [44| . 
The most massive galaxies are primarily ellipticals, while 
smaller galaxies are a mixture of ellipticals and spirals 
[45| . A more sophisticated Monte Carlo treatment would 
draw galactic bulge fraction from an observed distribu- 
tion as a function of galactic mass, and only correlate the 
SBH mass with this bulge component. Such a treatment 
is beyond the scope of this paper. 

This procedure for determining the SBH mass ratio 
q = 7712/7711 is considerably more elaborate than that 
used in the LZ papers [25|, l2g| . They assumed that the 
probability of having a SBH mass ratio between q and 
q + dq was simply F[q)dq, with the function F[q) given 
by Eq. ([2]) with c ~ 0.3 and rf = 1. This assumption 
seems to be based on a misinterpretation of Stewart et 
al. [33I, which defines F{q) as being proportional to the 
total number of mergers between a mass ratio of q and 
1. The LZ papers did not need to introduce an explicit 
cutoff at small mass ratios, because the total number of 
mergers remained finite under their assumption. This 
finite number of mergers allowed them to calculate the 
probability that the gravitational recoil would be in a 
particular velocity interval. Given the infinite rate of 
mergers predicted by Eq. ([1]) as m^ — >■ 0, we would pre- 
dict an infinite number of mergers with f = km/s in the 
absence of a cutoff, since v (x q^ as q ^>- 0. This cutoff is 
well justified physically because DM halos of arbitrarily 
small mass cannot host stars, let alone SBHs. 



B. Initial Spins 

The accretion and merger history of SBHs determines 
the magnitude and direction of their spins. In this section 
we discuss some astrophysical predictions on SBH spin 
magnitude and direction. 



1. Spin Magnitudes 

Here we present a short overview of the literature 
to show that SBH spin magnitudes are strongly depen- 
dent on the underlying assumptions about gas accretion. 
Because of these uncertainties, in this paper we sim- 
ply assume Xi = X2 = X for x ~ 0.5, 0.75, and 1.0, 
and we provide predictions for each of these three val- 
ues. Schnittman [20| showed that the amount of PN 
spin alignment is insensitive to the spin magnitude for 
X ^ 0.5 (see Fig. 11 in his paper). At large separations 
where |L| ^ \Xi\7 Eq. (3.5) of [2^ indicates that the con- 



straint for defining a spin-orbit resonance depends only 
on 0i, not Xi- We therefore expect that PN spin align- 
ment should not be particularly sensitive to our assump- 
tion that xi = X2i although additional numerical studies 
would be required to confirm this expectation. 

Volontcri et al. [ia] examined the distribution of SBH 
spin magnitudes using cosmological merger trees con- 
structed within the Press-Schechtcr formalism. They 
found that binary mergers alone lead to broad distribu- 
tions of the dimensionless spin magnitude peaked around 
X — 0.5, but that when accretion is included the spins 
reach much larger values. Accretion from geometrically 
thin disks [47[ leads to spin distributions sharply peaked 
around x — 0.998 [i^, as the Bardeen-Petterson effect 
14 1 aligns the SBH spins with the disks on a timescale 
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on which the spin magnitudes change. Here 770.1 is the 
radiative efficiency normalized by a typical value 0.1 and 
/Edd is the AGN luminosity in units of the Eddington 
value iEdd = 47rG'AfBHC777p/CTT (fTT-p is the proton mass, 
(Tt is the Thomson cross section for elastic scattering). 
Accretion from geometrically thick disks leads to much 
broader spin distributions, because spin alignment occurs 
on the longer Salpeter time. SBHs that accrete from a 
thick disk may also have a smaller maximum spin x — 
0.93, as ordered magnetic fields in the plunging region 
interior to the innermost stable orbit may extract angular 
momentum that would be advected by the SBH in the 
thin-disk case p9| . 

Volontcri et al. [4a] assumed that after each major 
merger, SBHs accreted continuously from disks with well 
defined angular momentum until their mass increased by 
an amount A?ti = 3.6 x IO^V^i^qMq, where 14,150 is the 
circular velocity of the host galaxy's DM halo in units 
of 150 km/s. Moderski and Sikora [501 proposed an al- 
ternative scenario, in which gas was accreted in small 
discrete episodes with random orientations with respect 
to the SBH spin. This "chaotic accretion" scenario leads 
to SBH spins that fluctuate about x = 0, and was origi- 
nally motivated by the conjecture of Wilson and Colbert 
51| that only radio-loud AGN (approximately 10% of all 
AGN) are powered by highly spinning SBHs. King et 
al. [52| suggested that the mass accreted in each of these 
discrete episodes might be determined by the mass of the 



accretion disk interior to the radius ' 



sg 



0.01-0.1 peat 



which the disk becomes self-gravitating. Star formation 
at r > Tgg would stir the gas flow, implying that each 
event would have essentially random orientations. Re- 
cent high-resolution N-body/SPH simulations of galactic 
accretion disks which include star formation support this 



suggestion [53|. These simulations follow gas flows from 
galactic scales of ~ 100 kpc all the way down to < 0.1 pc, 
and show that the orientation of the inner nuclear disk 
varies on '^ Myr timescales, that are comparable to the 
alignment timescale iai of Eq. ([5]) . 

Berti and Volonteri [5J| explored how chaotic accre- 
tion would affect SBH spin evolution in a cosmological 
context, also incorporating the results of NR simulations 
available at the time. The NR simulations implied that 
SBH spin distributions would now be broadly peaked 
about X ^ 0.7 in the absence of accretion, but that the 
full range of spin magnitudes < x < 1 was now possible 
once both standard and chaotic accretion were permit- 
ted. This differs from th e sp in magnitude distributions 
used in the LZ papers l23,[2a] , which as shown in Figure 7 
of [2g| are peaked around x — 0-8. These spin distribu- 
tions were determined from the N-body/SPH simulations 
of Dotti et al. [I^l, which were restricted to equal-mass 
binaries, and began with an initially uniform distribution 
of spin magnitudes Xi G [^A]- Although the < 10 Myr 
duration of these simulations is long compared to the 
alignment time iai on which the spin directions change, 
it is less than the Salpeter time iEdd on which the spin 
magnitudes change. One might therefore be concerned 
that these distributions have not converged to the value 
they should have at r^ ^ lOOOM, where the gravitational- 
wave driven stage of the inspiral begins. 



2. Spin Directions 

The residual misalignment of the SBH spins with their 
orbital angular momentum at r^ has been less thoroughly 
investigated than the spin magnitudes. We review some 
of these studies below. To encompass all possible scenar- 
ios with regard to astrophysical spin alignment prior to 
r,; = lOOOAf , in our simulations we will consider: 

1) Uniform symmetric distributions in cos^i and 
cos 6*2 drawn in the range Oi G [0°, 10°] (i = 1,2; 
strong alignment) and 9i £ [0°,30°] (weak align- 
ment). We will refer to these cases as the "10/10" 
and "30/30" scenarios. 

2) A uniform distribution in cos^i and cos 02 with 
asymmetric range for the angles: 6i G [0°,10°], 
02 e [0°,30°]. In this "10/30" scenario the primary 
is more strongly aligned with the orbital angular 
momentum. 

3) A uniform distribution in cos^i and cos 02 with 
6*1 G [0°,30°], 02 G [0°,10°]. This "30/10" scenario 
could result from the secondary orbiting close to 
the inner edge of a circumbinary disk, so that it 
ends up being more aligned than the primary with 
the orbital angular momentum. 

In order to justify these assumptions, here we review 
some literature on astrophysical predictions for the spin 
orientation resulting from SBH mergers. 



Dotti et al. [17[ performed simulations where one 
4 X IO^Mq SBH is initially at rest at the center of a 
lO^M© circumnuclear disk. The second SBH, also with 
7Ti,; = 4 X 1O^M0, began at an initial separation of 
r = 50 pc with orbital angular momentum that was al- 
ready aligned or anti-aligned with that of the circumnu- 
clear disk. This is consistent with the SPH simulations 
of Larwood and Papaloizou [55| , which suggest that the 
mass-quadrupole moment of the SBH binary will induce 
differential precession in an inclined circumbinary disk. 
Viscous dissipation will then cause this differentially pre- 
cessing disk to settle into the equatorial plane of the bi- 
nary. It is unclear whether this mechanism will operate 
at the large initial separation of the Dotti et al. sim- 
ulations before the SBH binary becomes gravitationally 
bound. Dynamical friction causes initially eccentric or 
retrograde orbits to circularize as the second SBH inspi- 
rals to a final separation r ~ 5 pc. Accretion during 
this inspiral aligns the spins of both SBHs to within 10° 
of the orbital angular momentum for a cold circumnu- 
clear disk, and to within 30° for a hot disk with greater 
pressure support, and hence more isotropic gas velocity 
dispersion [lOf. These simulations did not include star 
formation or cooling, and the internal energy of the gas 
particles was chosen to prevent gravitational fragmenta- 
tion. These choices suppress chaotic accretion and the 
dynamical instabilities found in Hopkins et al. |53] . that 
cause the direction of the angular momentum of the inner 
disk to fluctuate on < Myr timescales. These estimates 
should therefore be considered lower bounds on the resid- 
ual misalignment between the SBH spins and their orbital 
angular momentum. 

After the simulations of Dotti et al. [131 ®'^'^i but be- 
fore the gravitational-wave driven stage of the inspiral 
begins at r^ = lOOOA/, the SBHs open a gap and form a 
true circumbinary disk. High-resolution hydrodynamical 
simulations indicate that the gas surface density will be 
sharply truncated at radii less than twice the semi-major 
axis of the binary [5Q|. Mass can flow from the inner 
edg e of this circumbinary disk onto the individual SBHs 
p^ |57| . Accretion rates are generally higher for the less 
massive SBH, which is further from the center of mass 
and hence closer to the inner edge of the disk. This point 
is absolutely crucial for the PN spin alignment that is the 
subject of this paper, as it is the relative value of 0i and 
02 that determines whether gravitational recoils are sup- 
pressed or enhanced. One might expect that the higher 
accretion rate onto the smaller SBH would lead to closer 
alignment between its spin and that of the orbital angu- 
lar momentum. Alternatively, it is possible that the spin 
of the more massive SBH would be more influential in de- 
termining the direction of the warp in the circumbinary 
disk that torques the SBH spins. The relative value of 
01 and 02 in the presence of a circumbinary disk remains 
very much an open question. 

A flnal complication is the growth of the binary's 
orbital eccentricity through its interaction with the 
circumbinary disk [58|. Although dynamical friction 
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FIG. 3. Distributions of ejected remnants as functions of the binary mass ratio (top) and of the mass of the host galaxy 
(bottom). Left panels: v^j is computed using Eq. ^. Riglit: we use a constant ejection velocity Uoj = 1,000 km/s. In the 
asymmetric cases (10/30 and 30/10) thin lines refer to kicks computed using the initial spin distributions, while thick lines 
refer to kicks computed from the final distributions. 



damped the orbital eccentricity at the large separations 
of the Dotti et al. simulations [ITI , after the formation of 
a circumbinary disk the orbital eccentricity again grows 
to a limiting value Ccrit — 0.7 |59j . Gravitational radi- 
ation circularizes the binary after it decouples from the 
circumbinary disk [60[, but the binary could still have 



-0.73 



q 



considerable residual eccentricity glisa fx Af 
when the frequency of the fundamental GW harmonic 
reaches /lisa = lO"'^ Hz [5^,[6l|. We will neglect residual 
eccentricity during the inspiral and any effects it might 
have on PN spin alignment. 



IV. RESULTS 

We carried out a systematic study of the probabil- 
ity that gravitational recoils eject SBHs from their host 
galaxies. We generated a population of binaries with 
astrophysically motivated mass ratios using the Monte 
Carlo prescription described in Section Fill Al To quantify 
the effect of SBH spin magnitudes, we computed ejec- 
tion probabilities for three fixed values of xi = X2 = X'- 



X=0.5, 0.75, and 1. To account for all possibilities given 
the great uncertainty in the spin-alignment distributions, 
we considered the four cases described in Section [III B 21 
two "symmetric" cases (10/10, 30/30) and two "asym- 
metric" cases (10/30 and 30/10). 

Our results for x — 1 are summarized in Figure |31 
This figure has four panels: our simulated binaries are 
binned by mass ratio q in the top panels, and by the 
stellar mass M* of the larger host galaxy in the bottom 
panels. In the right panels we assume that a gravita- 
tional recoil greater than Woj = 1, 000 km/s is needed to 
eject a recoiling SBH from its host galaxy, while in the 
left panels we adopt a prescription that depends on the 
stellar mass Af*. The gravitational potential of a real 
galaxy is the sum of contributions from its stellar mass 
and DM halo; although the DM halo is more massive, 
the more concentrated stellar contribution is dominant 
for large elliptical galaxies [231 ■ Modeling elliptical galax- 
ies with the density profiles of Terzic and Graham J62l | , 
Gualandris and Merritt [63 find that escape velocities 
have typical values Uesc — 2.1(GA/*/i?e)^/^, where R^ is 
the effective radius that contains half the galaxy's total 



luminosity. Combining this with the observed relation 
Re « 1.2 kpc(M,/lOi°Af0)°o'^5 yigijs ^j^g approximate 

expression [63| 



l,154km/s 



Af, 



IQiiMp 



0.4625 



(7) 



This is the prescription we used to calculate escape ve- 
locities in the left panels of Figure [3] 

Each panel of Figure [3] has 7 histograms: the solid 
black curve corresponds to the total distribution of simu- 
lated binaries, while the 6 dashed and/or dotted colored 
curves correspond to the distributions of ejected SBHs 
given the 6 different spin-alignment distributions indi- 
cated in the keys to the figure. PN spin alignment can 
play an important role for asymmetric distributions, as 
can be seen by comparing the distributions of ejected bi- 
naries calculated from spins before (thin-line dashed and 
dash-dotted histograms) or after (thick-line histograms) 
the PN evolution from r, = lOOOM to r/ = lOM. As 
expected, the fraction of ejected binaries decreases in the 
10/30 case, while it increases in the 30/10 case. However, 
the kick reduction in the 10/30 case is more significant 
than the kick enhancement in the 30/10 case: in other 
words, PN resonances are more effective at reducing re- 
coils when the primary is more aligned with the orbital 
angular momentum than they are at increasing recoils 
when the secondary is more aligned with the orbital angu- 
lar momentum. For the symmetric distributions (10/10 
and 30/30), PN spin alignment has a marginal effect on 
the ejection probabilities For clarity (and to save com- 
putational time), in these symmetric cases we only plot 
the distribution of ejected binaries that results by ap- 
plying the recoil formula to the initial distribution at 
r^ = lOOOM. 

Comparing the left and right panels of Figure [3] allows 
us to understand the effects of the mass-dependent es- 
cape velocity of Eq. (O on the ejection probability. The 
most obvious effect is the steep decrease in the ejected 
fraction from near unity at M* < IO^'^Mq to less than 
10% for M* > 10"Mq. This is very natural: heavier 
galaxies should be more effective at retaining recoiling 
SBHs. Somewhat more surprising is that despite the con- 
stant escape velocity, the escape fraction also increases as 
M^, — >■ 1O^A^0 in the bottom right panel. This is a conse- 
quence of our decision to neglect SBHs with m2 < 10^ Mq 
that reside in galaxies with m* < 10^ Mq. Because of this 
choice, SBHs in galaxies with Af* ~ lO^Af0 only undergo 
comparable-mass (g ~ 1) mergers, elevating the fraction 
of large recoils since kicks are proportional to q^ /{1 + q)^ 
at leading order (see Eq. (11) of [26j). 

A second effect of the mass-dependent escape velocity 
of Eq. ([7]) is the possibility of ejecting recoiling SBHs for 



^ We have checked this statement by running PN evolutions in 
the 10/10 case. The "initial" and "final" distributions of ejected 
merger remnants are very similar. 



mass ratios as low as g ~ 0.1, as can be seen in the top left 
panel of Figure |3l The larger host galaxy can have a mass 
as small as IQ^^Mq for these mass ratios, corresponding 
to an escape velocity Wcsc — 400 km/s < 1,000 km/s 
by Eq. ([7]). This accounts for the nonzero escape prob- 
abilities for q < 0.1. We also note that in preparing 
these figures we have not evolved any PN inspirals for 
q < 0.1, even for those histograms labeled with N-l' that 
correspond to evolved distributions. This explains why 
the thin and thick lines corresponding to the initial and 
final histograms coincide for q < 0.1. The gravitational- 
wave inspiral time is inversely proportional to the bi- 
nary's symmetric mass ratio r] = q/(l -\- q)^ [60|, making 
it hard to preserve accuracy over long PN time evolutions 
for q < 0.1. However, Schnittman [20[ showed that spin- 
orbit resonances only become important inside a "reso- 
nance locking" radius nock/A/ w [(1 + q^)/(l - g^)]^. We 
are therefore justified in treating all binaries with g < 0.1 
(for which PN evolutions from r^ = lOOOAf to rj = lOM 
become problematic) as if the inspiral did not happen: 
that is, we can compute recoils of binaries with g < 0.1 
by simply applying the recoil formula at r = r^. 

We examine how the fraction /cj of ejected SBHs de- 
pends on the SBH spin magnitude x in Figured The two 
top panels in this figure were produced using the same 
distributions shown in the left panels of Figure [3] the 
fraction /cj is calculated by dividing the histograms for 
each spin distribution by the solid black histogram show- 
ing the total distribution. As expected, the ejected frac- 
tions are smaller for lower values of x, since the kick mag- 
nitudes are reduced while the escape velocities remain 
fixed. Ejections from host galaxies with M* > lO^^'^Afo 
(lO"A//0) are nearly eliminated for x ^ 0.75 (0.5). PN 
spin alignment remains effective for these lower spin mag- 
nitudes; for certain values of Af*, the ejected fraction 
is reduced by almost an order of magnitude. Careful 
scrutiny of the left panels of Figure 2] reveals that PN 
kick enhancement for the 30/10 case drives /ej above that 
of the 30/30 case, except at the highest values of M^, for 
the X = 0-5 and 0.75 distributions. Mergers involving 
the most massive host galaxies are dominated by small 
values of q, for which PN spin alignment is ineffective. 

It is also interesting to note that PN spin alignment is 
suppressed as g — )■ 1 (the thin and thick lines converge in 
this limit in the right panels of FigureH]) . Schnittman |2(j | 
showed that PN spin alignment increased as q increased 
from 1/9 up to 9/11 (see Figure 10 of Hol). Symmetry 
implies that PN spin alignment must vanish for equal- 
mass binaries, since the labeling of the SBHs is arbitrary 
for q = 1 and it is therefore impossible to distinguish 
6 1 < 62 from 9 1 > 62- This limiting behavior implies 
that the effects of spin alignment must be maximized for 
some mass ratio in the range 0.8 < g < 1, which is pre- 
cisely what we observe in the right panels of Figure SI 
These near-unity mass ratios are also expected to yield 
the largest kicks, which emphasizes the importance of in- 
cluding the effects of PN spin alignment in future studies 
of gravitational recoils. 
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FIG. 4. Fraction /oj of ejected remnants as functions of the mass of the host galaxy (left) and of the binary mass ratio (right). 
Here Vc] is computed using Eq. (O and we vary the dimensionless spin parameter X = Xi = X2- 



A final summary of our results is provided in Table U 
where we calculate the percentage of ejected binaries for 
each spin distribution. Although the precise amount of 
spin alignment depends on the value of Xi we see that 
kicks are suppressed (enhanced) by ~ 40% (20%) when 
the escape velocities are given by Eq. ([7]) . The kick sup- 
pression or enhancement is even more dramatic for the 
case of a constant escape velocity ^ej =1, 000 km/s. The 
ejected fraction is lower overall for this choice of escape 
velocity, which is larger than the mean value predicted by 
Eq. ([T])- The kick suppression or enhancement is greater 



because the effects of PN spin alignment are most pro- 
nounced at the high-velocity tail of the kick distribution, 
as shown in the left panel of Figure [5] 



V. DISCUSSION 

We demonstrated in previous work [22| that the align- 
ment of SBH spins with each other during the PN stage 
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TABLE I. Percentage of ejected binaries for different values 
of the spin and different combinations of the angles {9\, 62) ■ 
A column header such as "(10/30)^" means that the proba- 
bility was computed considering the spin distribution before 
PN evolutions, while "(10/30)/" means that the probability 
wEis computed applying the recoil formula at the end of the 
PN evolution. In the top rows we assume that the escape 
velocity Vej is given by Eq. (0 ; in the bottom rows we assume 
a constant v^ = 1,000 km/s. 



of the inspiral can dramatically suppress the predicted 
distribution of gravitational recoils. This mechanism can 
reinforce the kick suppression that results from the align- 
ment of the SBH spins with their orbital angular momen- 
tum through interaction with a circumbinary disk at an 
earlier stage of the inspiral [ia| . Recent NR simulations 
|25l [26j showed that gravitational recoils are maximized 
not in the previously claimed "superkick" configuration 
in which the spins lie in the orbital plane, but in a new 
"hang-up" configuration in which the angle between the 
spins and orbital angular momentum is ~ 50°. The pri- 
mary conclusion of this paper is that kick suppression due 
to PN spin alignment remains highly effective despite the 
new dependence of the gravitational recoil on SBH spins. 

It is difficult to make very robust quantitative predic- 
tions about the magnitude of this effect, as the recoil 
distribution is highly sensitive to SBH spin distributions 
that are theoretically uncertain and poorly constrained 
observationally. Some SPH simulations have shown that 
SBH spins can grow to large magnitudes and rapidly align 
with the orbital angular momentum through coherent ac- 
cretion from a massive circumnuclear disk [l9[ , but other 
simulations that include gas cooling and star formation 
indicate that the direction of the angular momentum of 
the accreted gas will vary chaotically on timescales com- 
parable to the alignment time, leading to smaller spin 
magnitudes and lar ger misalignments with the orbital 
angular momentum [53|. Our results show that the pre- 
dicted fraction of recoiling SBHs that are ejected from 
their host galaxies can vary from < 10~^ to ~ 1/3 de- 
pending on the adopted distribution of SBH spins. Al- 
though somewhat frustrating from a theoretical perspec- 
tive, the strong dependence of the ejected fraction on 



SBH spins implies that observational studies of recoils 
may place promising constraints on the highly elusive 
SBH spin distribution. 

Our work in this paper reveals several of the issues that 
must be addressed before we can predict recoil distribu- 
tions with confidence: 

1) We need to determine how effectively realistic cir- 
cumbinary disks can align SBH spins with their or- 
bital angular momentum. Such circumbinary disks 
may not remain geometrically thin, and may grav- 
itationally fragment or collapse into stars. These 
possibilities could have significant quantitative ef- 
fects on the SBH spin distribution at the onset of 
the PN stage of the inspiral. 

2) We need to develop a better understanding of the 
interaction between unequal-mass SBHs and a cir- 
cumbinary disk. One could argue that the spin 
of either the primary or the secondary might align 
more efficiently with the orbital angular momen- 
tum. Whether PN spin alignment leads to kick 
suppression or enhancement depends crucially on 
the relative values of 9i and 62 ■ 

3) We need to carefully consider how the merger rate, 
mass-ratio distribution, and escape velocity depend 
on host-galaxy mass. Our results shown in Figs. [3] 
and m demonstrate the sensitivity of the ejected 
fraction to assumptions about these galactic prop- 
erties. 

Although challenging, steady theoretical and observa- 
tional progress is being made on all of these issues. Once 
astrophysicists can provide more accurate predictions of 
SBH spin distributions at r^ - lOOOAf , when SBHs de- 
couple from their circumbinary disks, relativists will be 
able to evolve the binaries through merger and more ac- 
curately predict gravitational recoils. Including the PN 
stage of the inspiral will be an important ingredient in 
this grand theoretical undertaking. 
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